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The Nambu-Jona-Lasinio (NJL) model is one of the most frequently used four-fermion models in 
the study of dynamical symmetry breaking. In particular, the NJL model is convenient for that 
analysis at finite temperature, chemical potential and size effects, as has been explored in the last 
decade. With this motivation, we investigate the finite-size effects on the phase structure of the 
NJL model in D = 3 Euclidean dimensions, in the situations that one, two and three dimensions are 
compactified. In this context, we employ the zeta-function and compactification methods to calcu- 
late the effective potential and gap equation. The critical lines that separate trivial and non-trivial 
fermion mass phases in a second order transition are obtained. We also analyze the system at finite 
temperature, considering the inverse of temperature as the size of one of the compactified dimensions. 

^ ; PACS number(s): fl.IO.Kk; fl.IO.Wx; If.30.Rd 
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The last decades witnessed significant investigations on the phase structure of quantum field theories, in particular 
on the chiral symmetry phase transitions in Quantum Cromodynamics (QCD). However, due to its complex structure, 
effective models have been largely employed to simplify that analysis. Among them, one of the most frequently used 
is the four-fermion theory known as Nambu-Jona-Lasinio (NJL) model 0. The NJL model is specially convenient for 
the investigation of dynamical symmetries when the system is under certain conditions, like finite temperature, finite 
chemical potential, external gauge field, gravitation field and others 

Finite-size effects on the phase transitions of four-fermion models have also attracted a great interest 0, 0, Q ■ This 
question emerges when the system has a finite size and it is not clear if it is large enough to apply the thermodynamic 
' limit in a usual way; frequently it is necessary to take into account the fluctuations due to finite-size effects. In 
particular, Ref. performed a numerical investigation of a three-dimensional four-fermion model in a finite-size 
' scaling analysis, where the finite-size effects act as an external field. On the other hand, Ref. studied the NJL 
. model in the framework of the multiple reflection expansion, where in terms of a modifled density of states finite-size 
Q_i' effects are included (see also Ref. 1^). The critical temperature was suggested to decrease as the system is reduced. 
D In this paper, we investigate finite-size effects on the dynamical symmetry breaking in a different way. We study 

the Euclidean three-dimensional NJL model in the framework of zeta-function and compactification methods _9J . This 
^ I procedure in principle allows us to explore the mentioned model with one, two or three compactified dimensions with 
: '"j • antiperiodic boundary conditions ^3 and compare their effects in the phase diagram of the model. With the choice 
r^J of all dimensions being spatial, the system is considered confined between two parallel planes a distance L apart, 
^ 1 confined to a infinity cylinder having a square transversal section of area L^, and to a cubic box of volume L^, for 
■ ■ ■ ' one, two and three compactified coordinates, respectively. At finite temperature, we associate one of the compactified 
coordinates to the range [0, /3], where /3 is the inverse of the temperature T. In this setting, we calculate and determine 
analytically the size-dependence of the effective potential and gap equation. Phase diagrams, where the symmetric 
and broken phases are separated by a size-dependent critical line, are obtained. 

This paper is organized as follows. In Section II, we briefly review the NJL model and obtain the effective potential 
in mean-field approximation. In Section III, we apply the zeta-function method to derive the effective potential. After, 
Section IV is devoted to analyze the size-dependent gap equation. The phase diagrams are shown and discussed in 
Section V. Finally, in Section VI are given a summary and concluding remarks. 
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II. THE FORMALISM 



NJL model as an effective theory of QCD 



QCD is the theory of strongly interacting matter, which is constituted of quarks and gluons. In the study of 
dynamical symmetry breaking, the attractive interaction between quarks and antiquarks, coming from complicated 
processes of exchanges of gluons, are often replaced by an effective interaction between them. This interaction 
generates the quark-antiquark condensation in the vacuum when the interaction exceeds a critical strength, and 
this condensation is responsible for the dynamical symmetry breaking. Thus, theories incorporating this mentioned 
symmetry, such as the NJL model, are very useful in this analysis. 

The NJL model can be obtained heuristically as follows. Inspired by the QCD, let us consider the following 
approximation for the action which describes the interactions between quarks mediated by exchange of one gluon. 



A. 



exchange 



(1) 



where the spinors q and q represent the quark and antiquark fields carrying Nf flavors and Nc colors, A^ is the gluon 
field, A° (a = l...(A^^ — 1)) are the generators of SU(7Vc) and M is the quark mass matrix. The quantity D^^ix — y) 
is interpreted as the effective gluon propagator. Note that there is no self-gluon interaction term. The functional 
integration over the gluon field in Eq. leads to 



A. 



exchange 



= d 



q (i^-m) q + ^ d'y {qX^^^q) {x)DZ{x - y) (qXl^.q) (y) 



(2) 



Next, we assume that the interaction obeys to the Feynman gauge, and moreover it is a contact interaction, which 
means D^^^ {x — y) oc GSabg'^^ 5{x — y) . Furthermore, we can use the Fierz identities for color, flavor and Dirac spaces 
0, El , to rewrite the interaction part of the Lagrangian in the action explicited in Eq. Q as (considering only the 
qq-current). 



A} ® K^, with A} = [l/y/N}, xf^ ; a 



Cm (X G {qk'^q) (qA^q) , (3) 
— 1, ■■■jN'f — 1 being the generators of U{Nf), and = 



where A" = Ic 

Eq. 10) characterizes the interaction part of the Lagrangian of the NJL model. It simplifies a gluon exchange 
interaction between quarks and antiquarks to a local interaction, i.e. a four-quark point interaction. 



B. The effective potential 

In the previous subsection we have heuristically derived the NJL model from a simplified version of QCD. To analyze 
With its dynamical symmetry phase diagram, we will begin by deriving its effective potential. For generality, we will 
consider the NJL model as an independent model by itself: the coupling constant is considered to be and arbitrary 
parameter, without any a priori relation with QCD aspects. Also, the dimension of the system is considered here 
arbitrary (2 < Z? < 4), as well as we limit ourselves to mesonic scalar and pseudoscalar modes. 

We will treat with the colorless and massless NJL model, 



C = q i ^ q 



E 

1=0 



(9*75 A* g) 



(4) 



where the spinors q and q here carry N flavors, and the matrices A' are the generators of the U{N), A' = 

(//\/]V,A") ; a = l,...,iV2~l. 

In order to study the phase structure of this model, it is convenient to perform the bosonization. We introduce 
auxiliary fields M* and 0' through the following way. 



VA'P S{M' - qVq) 



VAPVcj)' exp 



i / d'^xiM' - qr'q)<j)' 



(5) 
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where (j)^ = (cr*, tt') and = (A*, 175 A*). In this sense, cr* and tt* are fields corresponding to the scalar and pseudoscalar 
bilinear forms, respectively. The use of auxiliary fields in the generating functional of the model introduced in Eq. 
(jH) yields 



^ cx / T)qT)qT>(j)' exp{i / d^x 



2G 

where the source terms have been omitted. Thus, the integration over the Grassmann variables q and q results in 



(6) 



j P?i'exp{iA//}, (7) 

where 

is the effective action associated to the model described in Eq. (0} with one-loop corrections; is a constant with 
dimension of mass. Here Tr means the trace over coordinate, Dirac and internal spaces. 

In the above context we restrict to be independent of the spatial coordinates; moreover, we impose tt* = 0, 
ct" = (a = 1, ...,iV^ — 1); only the scalar part of cf) which is diagonal in flavor space, cr°, can assume non-vanishing 
values. Thus, it is convenient to define A with a assuming the role of dynamical fcrmion mass, its non-trivial 

value meaning that the system would be in the broken phase. 

Taking into account the above considerations and performing the Wick rotation in the x'^-coordinate, we can use 
Eq. © to obtain the the effective potential, 

lc,„ = -^ = ^ + a,M, (9) 

where V is the volume, and 

with h being the dimension of the Dirac representation. 

As our interest is to describe the finite-size effects on the phase structure of this model, it is convenient to use the 
zeta- function regularization techniques, described below. 

III. THE ZETA-FUNCTION REGULARIZATION APPROACH 
A. The zeta-function 

This method is based on the fact that one can define a generalized zeta-function from the eigenvalues a; of a 
differential operator A of order 2 on a d-dimensional compact manifold (a bounded manifold) , 

a(s) = E"r, (11) 

i 

valid for Re s > ^. For the other values of s, we must perform an analytical continuation (for a review, see [9I 
In particular. 



ds 



= -^Ina^ = -TrlnA. (12) 



Therefore one can associate the operator A to the one appearing in the argument of the logarithm in Eq. I|10|l . 
allowing to rewrite Ui{a) as 

Ui{a) = ^[CA{0) + ^nti\A{0)]. (13) 
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B. Compactification 

At this point we will take into account the finite size effects. Considering the system restricted to have d < D 
compactified dimensions, we can rewrite the coordinate vectors as xe — {y, z), where 

y = (yi = x^, ...,?/„ EE x'^) , z = {zi = x'^+'^,...,Zd = Xe) , (14) 

with the Zi components defined in the domain Zi £ [0, Li] and n = D — d. The compactification of the coordinates 
makes the fc^-components of momenta to assume discrete values, 

^2 ^ = 0,1,2,...; (15) 

where Ci = ^ {i — 1,2, ...,c?) for antiperiodic boundary conditions. Notice that we may associate a given Li to the 
inverse of temperature (3 = 1/T, i.e. Li = f3; in this situation the chemical potential fiQ also could be introduced 
through the rule Ci = i - 

The prescription given in Eq. H15|) makes the zeta-function in (|12|) be given by 

ni , . . . ,71^ — — oo ^ ^ 

where Vn is the n-dimensional volume. The techniques of dimensional regularization can be used to perform the 
integration over the fcy-components, yielding 

U is; {aj , {cj) = ^^^^Yf {s ~ |; {aj , {cj) , (17) 

where (y\ {ai\ , {c^}) is the generahzed Epstcin-zcta function, defined by 

Yf (y-Mi\ M}) = E 

Til ... . ,nd— — oo 

being = notice that Y^ is valid for Re > ^. 



ai (ni + cxf ^ h ad (rid + c^)^ + i 



(18) 



C. The pole structure of the zeta-function 

According to Eq. (|13() . in order to obtain the effective potential, we must perform the derivative of the 
(s; {ci} , {ci}) at s = 0. However, we note from Eq. H17|) that there is a different pole structure coming from 

\i ? ' at s = for n even or not. So, we are forced to consider these different situations. Furthermore, the function 
Y^^ also has a different behavior for d even or odd, as it is discussed in Ref. [l3|. In the case of d odd, Y^^ has poles 
at 1/ = |, I — 1, i and _ (^!±ll (/ g Nq), while for d even they are at z/ = |, | — 1, 1. Nevertheless, as it will be 
shown in the next section, for the study of the phase diagram we take 0; in this limit, Y^ has pole only at |. 

2 

This fact assures that for further calculations the pole structure of Y^ is not relevant. 
Then, from Eq. H17|l we have to deal with the following situation. 



C:i(0; {aj,{cj) = 



(47r)"/2 s^o 

r(^-t) d 
r (s) ds 



T{s) 

(Yf (.-|;{aj,{c,} 



(19) 



After the necessary manipulations, we obtain 



C;i(0;{aj,{c,}) = 



21\ 
2 • 



s- ■^■Aai},{c^} 



(47r)"/2 

yf [7 + V'(f + l)]}, 



for n even, 



(20) 
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where 7 is the Euler-Mascheroni constant and tp (k) is the digamma function, and 

Vn - 



CAiO;{a^}Ac^}) 



(47r)"/2 



Yf (-^;{aj,{c.}), 



for n odd. 



(21) 



Due to the discussion of the pole structure above, notice that Eq. 12111 is vahd for d even or in the hmit ct^ ^ 0. 
Now, coming back to Eq. p3|l . we have two possibihties for Ui{a): the use of Eq. H2U|) gives 



Ui (ct; {a,,} , {c^}) = 



h (-1) 



for n even, 



while with the help of Eq. (|21|l we obtain 



for n odd, 



(22) 



(23) 



valid for d even or in the limit ct^ ^ q. 

Now we need to perform the analytical continuation of the generalized Epstein zeta-function YJ {O'i} i {ci})- 

Making use of Mellin transforms and relations between Jacobi theta-functions ^ , we find 



(:.;{aj,{cj) = 



7r2 



d 2 



E 



1 = 1 J 



X cos{2TmiCi) . . .cos{2TmdCd)K^_d \2ti(t 



(24) 



where Z means the set of non- vanishing integers. 

For small values of cr^, cr^ << 1, it is possible to use the binomial expansion for Y^ defined in Eq. p8|) to expand 
it in powers of the cr-field , obtaining the expression 



Yf {q; {a J , {cj) ^ -^nq,m (q + j; {«J , {cj) a'^, 
3=0 ^' 



(25) 



where T(q,j 



- p(^) , and 



Yd {ly; {ai} , {ci}) ^ ^ ^ai{7ii + cif + ... + ad{nd + Cdf 
is the homogeneous generalized Epstein zeta-function, with Zq being the set of integers. 



(26) 



IV. THE GAP EQUATION 
A. Case without boundaries 



In this section we focus our attention to the gap equation, which will allow us to analyze the phase structure of the 
model. It is obtained by minimizing the effective potential with respect to a, 



d 

g^Ueff (cr;{a,},{cj) 



= 0, 



(27) 



where m is the dynamical mass of the fcrmion. 
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For completeness, we briefly describe the case without boundaries (for reviews see Refs. [Ullllll). This situation 
means d — 0, and therefore n — D. So, it is more useful to use the formula H23I) of Ui for n odd, 



Ui (cr; d = 0) 



D{iTT)D/2 



r 1 



D 



n 



(28) 



Notice that Ui diverges for D even [l^. In particular, for D — 2, four-fermion theories are renormalizable, and the 
renormalization is implemented by imposing the following renormalization condition. 



da 



j^Ueff (a;d = 0) 



(T—fl 



1 

Gr ' 



which, with the help of Eqs. © and yields 

1 1 h{D-l] 



r 1 



D 



,D-2 



Gr G (47r)^/2 ■ V' 2 
Thus, we can write the renormalized effective potential, valid for 2 < _D < 4, as 

Hence, the non-trivial solution of the gap equation is given by 

where we have defined the dimensionless coupling constant Gi = h^^^Gr, and 

(4^)^/2 



D-2 



Go 



hii-D)r{i-§) 



(29) 



(30) 



(31) 



(32) 



(33) 



From Eq. (|32|l . we can infer that for Gi < Go the allowed value of the mass m is trivial, while for Gi > Go the mass 
acquires a nonzero positive value. In this context, the constant Go acts as a critical value which Gi must exceed to 
have dynamically generated fermion mass. 



B. The presence of boundaries 



In the case of presence of boundaries, we can use the same systematics as above, and therefore determine the 
finite-size contributions to the expression for Gi, with the only difference that we must analyze the gap equation in 
Eq. (|32l) with its last term in right hand side in a different way. 



1 



1 



d 



n - n - D-2 TT^i (""'•>^°»J'''!^'^'J') 
Gi Go m^^ ^ oa 



(34) 



where m is the fermion mass with the system in the presence of boundaries. 

We can calculate the derivative of J7i with respect to a directly from Eq. (|13() . by acting it on the zeta-functions. 
Taking into account the different possibilities for [/i, then the gap equation is written in the different situations, 



1 

'g[ 



1 



2-D 



(-1) 



>- rl 



Go Ki(47r)«/2 (I _ 1) 



+ i;{a,},{cj 



HE) 



for n even; 



while 



2-D 



Gi Go yrf(47r)"/2 



(l - f ) Yf (-f + 1; W^} , , for n odd; 



(35) 



(36) 
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(remembering that Eq. 136() is valid for d even or in the hmit cr^ — > 0) and finally we have obtained, 

1 1 hfi^-^ 



Gi Go Vd 



FPy7-(l;{aJ,{cJ), forn^O, (37) 



in which we have considered only the relevant terms for further calculations; FP means the finite part of YJ" . 

Now, we are able to determine the critical value of the coupling constant Gi with the corrections due to the presence 
of boundaries, just by taking the fermion mass approaching to zero in the gap equation. In this context, we can use the 
expansion in Eq. H25|l for m ^ 0; taking only the first term of this expansion, then FJ" reduces to the homogeneous 
generalized Epstein zeta- function Yd- However, as it was remarked in the previous section , the pole structure of Yd 
is simpler than that of YJ" ; therefore Eq. H36() becomes valid for general d. 

We can construct analytical continuations for Yd to write it in terms of Bessel and Riemann zeta functions, by 
considering a generalization of recurrence formulas in jSj, 

Yd {v, {aj , {c,}) = ^S^^^Yd-i [v-]^- {a,^d} , {c,^d}) + ^W^d (^^ - \; {«J , {cj) , (38) 
where the set {aj^d} means that the parameter ad is excluded from it, and 

Wd im {a^} , {c.}) = ^ E '^os i^TTUdCd) ( ' y ) ' (^Xd-i) , (39) 



with Xd-i = Y X]fe=i '^k {nk + Cfe)^ and K^{z) being the modified Bessel function of second kind. 

Notice from Eq. 1)3811 that whatever the sum one chooses to perform firstly, when Ci — Cj the Li <-> Lj symmetry is 
lost |l5j |: in order to preserve this symmetry, we adopt here a symmetrized summation generalizing the prescription 
explicit in for = 0. 

In the next section we will consider particular cases. 

V. PHASE BOUNDARY 

In this section we will analyze the expressions (|36|l - (|37|) for m ^ 0, which defines a critical surface L^-dependent 
for the phase diagram of fermion mass. In the first three subsections, we consider the situations where all three 
dimensions are spatial coordinates. In the last subsection, we analyze the system at finite temperature, associating 
one of compactified dimension to the range of the inverse of temperature. 

A. Case d = 1, n even 

We now consider the simplest case of the compactification, namely when only one coordinate is compactified. This 
situation was in part studied in Ref. |l2j considering Li = /3; for completeness we will describe it here. 
For d = 1, Eq. is written as 

Fi {i^; ai = a, ci = c) = a'"' [C{2v, c) + C(2i^, 1 - c)] , (40) 

where C(s, c) = X^fc^o + * is the Hurwitz zeta-function. The use of this last equation in (|35|l for ?fi ^ yields 

Gi Go L(47r)— ("^j' 

[C{D~3,c) + aD-3,l-c)]\. (41) 



D-3 

a 2 



In/i^ - ^ ( ^ ^ ^ ) - 7 - Ina 



Employing the expansion formula for the Hurwitz zeta-function. 



Inr(c) - iln27r 



(42) 
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Eq. glj becomes, for D = 3, 



11 
Gi Go 

Notice that for c = ^ Eq. H43|l reduces to 



■In 



27r 



^TL^l Vr(c)r(i-c) 



(43) 



1 

G~i 



1 



(44) 



where Ai = = 0.221. 

Noting that using H32|l in H44|l . we obtain the formula Lm = 21n2. It can be seen that the replacement of the length 
L by the inverse of temperature (3 gives the equation Tc — 2hi2 ^ which can be identified as the critical temperature 
of a second order phase transition, as has been remarked in Ref. [||. 



Furthermore, yet with L = f3 and c 



_ 1 _ 

2 
1 



1 

Go 



where is the chemical potential, Eq. (|43|l becomes 



■In 



2 cosh 



(45) 



The use of Eq. (|32|) in H45|) generates the following formula for the critical temperature and the chemical potential. 



13 m 



In 



2 cosh ^ 



, which agrees with that obtained in Ref. (ij 



B. Case d = 2, n odd 

We now focus our interest in the case of two compactified coordinates. Taking d = 2 in Eq. (|36f) with m ^ 
generates 



1 1 



Gi Go V"2(47r)^ 



Y, - 



;ai,a2,Ci,C2 



(46) 



The analytical continuation of this equation can be performed by using H38() in its symmetrized version, 140() and 142|l , 
allowing us to rewrite Eq. 14()|l for _D = 3 as 



1 1 



27r 



Gi Go Ml^27r|2^ Vr(c,)r(l-c,) 

with j ^ i. 

Considering q = i, we obtain [j ^ i) 



oo oo 



i—l rii—l nj— — oo 



(47) 



1 



1 1/1 



Gi Go TT I \2LifL 2L2H ^ 
In particular, for Li = L2 = L, we have 

where 



Gi " Go L^i' 



00 CXD 



77- 77 r J ( J \ 



rii — 1 n2 — — 00 



"2 



9^ 0.257 . 



(48) 



(49) 



(50) 
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C. Case d — 3, n = 

We now pay our attention to the case where three dimensions are compactifieci. Taking d = 3 in Eq. H37|l. we 
obtain for m = 0, 



1 _ 1 2 



FP Ys (1; ai, a2, as, ci, C2, C3) . 



(51) 



The analytical structure of the function Y3 (1; ai, 02, 03, ci, C2, C3) can be obtained from the general symmetrized 
recurrence relation given by Eq. H38|) ; explicitly, one has 



Ys (l;ai, 02,03, ci,C2,C3) = - ^ 



(l + eyfc) 1 



i.j,k—l 



1 



An 



Y2 -;aj,ak,Cj,Ck + —W3 -; oi, 02, 03, ci, C2, C3 , (52) 



where Sijk is the totally antisymmetric symbol. With the same expression for Y2 as used in Eq. H46|l . we can rewrite 
(jFT|l as 



1 (l+£^.. 



1 



Gi Go inV^n 



L 



27r 



2 vr(cfe)r(i-cfe) 



00 00 



Lj cos {2TTnjCj)KQ (2TrLjnjTk) + {Lj <-!■ Lfe) 



rij— 1 nfc— — 00 
3 



^ {l+Eijk) Li ^ cos(27rniCi) ^ 

i,j,k=l ni = l 



rij ,nfc— — 00 



TliLi 



T 



J k 



Ki {2TrniLiTj k) 



(53) 



where T,, 



— — ^ — I?. — 



If we restrict ourselves to the situation where Ci = ^, then we obtain the following equation that describes the 
critical coupling. 



^ . V ii±£^Lj^ln2 

Gi Go SnVafi | . ^ ^ 2 ^ 2 

00 00 

- L,J2 (-l)"'^o {2^L,n,n) + (i, ^ Lfc) 



TLj—l Ilk — — 00 



J2 (i+£,,.)i.E(-i)"' E 

i,j,k=l ni = l 



TliLi 



/-iTi I 2TTniLiTj k 



(54) 



where Ty... = ^1. 

In particular, for Li = L2 = L^ = L, we get 



where 



Gi Go L^,' 



(55) 



A3 = i^-^f: £ (-l)-i^o(2.ni 

m — l 712— —00 



^^2 



ni 



"2,"3 = — 00 



i=2 V"* 2; 



\ 



fii H 

2 



= 0.278. 



(56) 
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In Fig. ^is plotted the critical coupling constant Gi as a function of x = {Lfj,)~^ representing Eqs. (|44|l . H50|) 
and H55|l. These contexts correspond, respectively, to the system in the following scenarios: confined between two 
parallel planes a distance L apart; confined to an infinitely long cylinder having a square transversal section of area 
L^; and to a cubic box of volume L^. We can see that as L decreases, Gi increases monotonically, meaning that 
to keep the non-trivial mass phase with compactified coordinates, it is necessary a stronger interaction than that in 
non-compactified space. 

In addition, from Fig. Qit is easy to see that there is a critical value Lc below which there exists no more non- 
vanishing fermion mass, since the coupling constant tends to infinity at this point. We have the following values for 
Xc = {Lcfi)~^ = 1.44 for d — 1, 1.24 for d = 2 and 1.14 for d = 3. Thus, the cubic box gives the greatest value for 
Lcj which means that the increasing of the number of compactified dimensions forces a stronger interaction to stay 
in the region of non-trivial mass. 



14 
12 
10 
" 8 
6 
4 



0.2 0.4 0.6 0.8 

X 

FIG. 1: The critical coupling constant Gi as a function of x = {Lfi)~^ for the three situations: dashed, dotted and solid lines 
represent Eqs. (1441 . II5U1 and (1551 . which are the cases of d = 1, 2, 3 compactified dimensions, respectively, L being the length 
in each case. For each situation the non-trivial mass phase corresponds to the region above the corresponding line. 




D. System at finite temperature 



Now we turn our attention to the system with compactified dimensions and at finite-temperature. First, let us look 
at the situation in Eq. (|48|l . which has d = 2. Taking Li = L and L2 = f3, we obtain the expression for Gi which 
is dependent of the temperature effects and of the size of the compactified dimension. However, we also can use Eq. 
obtaining 



m — 1 n2 — — 00 



^ 00 00 ^ J 

ni = ln2=— 00 ^ 



^2 



(57) 



The corresponding phase diagram in the x — T plane of Eq. (|57|l is plotted in Fig. |21 where x = {Lm)~^ and 
T — {fim)^^ . Notice that the critical line separates the non- vanishing fermion mass phase below the line and the 
trivial fermion mass phase above the line. As we can see, the size decreasing reduces the critical temperature. Also, 
the diagram has a manifest x <-^- T symmetry. 

Furthermore, in the L ^ 00 limit, which in the diagram means x ^ 0, the variable T reaches the usual critical value 
Tc = 1/(2 In 2). Analytically, this results from Eq. H57I) . in which only the first term with Bessel function survives, 
and the summation over 712 is replaced by an integration, yielding 



I tt (-!)-«. (2,f,„ 



1 

n2 -\ — 
2 



L^oo 4 



Jo 27r 

.i-hm(l-2-)C(.s) = -^ln2. 



(58) 



11 



FIG. 2; The phase diagram corresponding to Eq. H57|l (finite T and one compactified spatial dimension) in the x — T plane, 
where x = {Lm)~^ and T = (/3m)~^. The non-zero mass phase corresponds to the region below the curve. 

Hence, the replacement of Eq. H58|) in H57I) in the hmit L — > c» reproduces the correct critical value 
Tc = 1/(2 In 2) = 0.72. Moreover, since the diagram has a x ^ T symmetry, the same argument is apphed 
to L in the limit f3 ~* oo. 

Now we analyze the case with three compactified dimensions, by using Eq. (|32|l in Eq. H54|) and taking Li — L2 ^ L 
and L3 = 13. Then, we obtain the following equation for critical line. 



TO 



L + ^)ln2-— y y (-l)"ii^o(2^^„i 

^ r I r / r — 1„2= — OO ^ 



n2 



. 00 00 / J 



3/3^ 
8 



m — 1 n2— — 00 

00 



n2 



J2 E i~^rKo(2nn, 



m — 1 712— —00 



"2 



"2, "3=— 00 



LI3 . 



m — 1 



(59) 



The corresponding phase diagram of Eq. (|59|1 in the x — T plane is plotted in Fig. O As we can see, there is no 
X ^ T symmetry in the case of the critical line of Eq. 159|) . Besides, the critical temperature has a faster decreasing 
with the size reduction than the case of Fig. |2| 

It should be noted that, by employing the same systematics used to get Eq. the x — > gives again Tc — 

1/(2 In 2) = 0.72. On the other hand, the T ^ limit generates x^ = 0.62. 



VI. CONCLUDING REMARKS 



In this paper, we have discussed the finite-size effects on the phase structure of the Euclidean three-dimensional 
NJL model in the frame of zeta function and compactification methods. 

In the leading order of the -^-expansion, we analytically derived the expressions for the renormalized effective 
potential and gap equation for different situations of compactified coordinates without any further approximation. 

Taking first the situation in which the dimensions are spatial and with antiperiodic boundary conditions, we 
estimated the strength of the critical coupling constant for the cases of one, two and three compactified dimensions 
with the same length L and compare them. We found that increasing of the number of compactified dimensions 
requires the increasing of the strength of the interaction to remain in the region of non-trivial mass. 
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FIG. 3: The phase diagram corresponding to Eq. H59^ (finite T and two compactified spatial dimensions) in the x — T plane; 
X — {Lm)~^ and T — {f3m)~^ . The non-zero mass phase corresponds to the region below the curve. 

We also investigated the system at finite temperature, which formally is done by associating the size of one of 
compactified dimensions to /3 = 1/T. We have obtained the phase diagram in terms of temperature and inverse of 
the size and noted that the phase transition is easier (bigger values of L) for the situation with the largest number of 
compactified coordinates. 

It is worthy to mention that we have obtained general expressions for the effective potential and gap equations 
(see Eqs. (I23ll - (|22ll and H36|l - I|37() '). for different numbers of total and compactified dimensions and also for different 
choices of the c^-parameters that appear in Eq. H15|) . This approach in principle allows us to study other cases, as for 
example the phase diagram at finite chemical potential as well as non-trivial topologies. 
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